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We report a measurement of the B — > iriv branching fraction based on 211 fb _1 of data collected 
with the BABAR detector. We use samples of B° and B + mesons tagged by a second B meson recon- 
structed in a semileptonic or hadronic decay, and combine the results assuming isospin symmetry 
to obtain B(B° -> -k~£ + v) = (1.33 ± 0.17 8ta t ± 0.11 syst ) X 10~ 4 . We determine the magnitude of 
the Cabibbo-Kobayashi-Maskawa matrix element \V u t \ by combining the partial branching fractions 
measured in ranges of the momentum transfer squared and theoretical calculations of the form fac- 
tor. Using a recent lattice QCD calculation, we find \V u b\ = (4.5 ± 0.5 s t a t i 0.3 sys t^Q^5PF ) x 10 _J , 
where the last error is due to the normalization of the form factor. 

PACS numbers: 13.20.He, 12.15.Hh, 12.38.Qk, 14.40.Nd 



The magnitude of the Cabibbo-Kobayashi-Maskawa 
matrix [l[ element V u b is a critical constraint on the Uni- 
tarity Triangle. Our knowledge of \V u b\ comes from mea- 
surements of the b — > ulv decay rate, where the hadronic 
system in the final state can be reconstructed either in- 
clusively or exclusively. The precisions are limited by the 
uncertainties in the non-perturbative QCD calculations 
that are used to extract IK^I from the measured decay 
rates. It is therefore crucial to pursue both the inclusive 
and exclusive approaches, which rely on different theo- 
retical methods, and to test their consistency. 

The rate of the exclusive decay B — ► ir£v (£ — e or \i) 
is related to jV^I through the form factor f+(q 2 ), where 
q 2 is the momentum transfer squared. Measurements 
of the B — > tt£v branching fraction have been reported 
by CLEO 0, BABAR H, and Belle Q. In this Letter, 
we report a measurement in which B — > tt£v decays are 
searched for in T(4S) — * BB events that are identified by 
reconstruction of the second B meson (Stag)- The tech- 
nique, which was also used in allows us to constrain 
the kinematics, reduce the combinatorics, and determine 
the charge of the signal B. The result is an improved sig- 
nal purity at the expense of the efficiency compared with 
the traditional measurements in which only the signal 
B meson is reconstructed. We perform two analyses in 
which -Btag is reconstructed in semileptonic and hadronic 
decays, respectively, and combine the measured partial 
branching fractions AB in ranges of q 2 with the recent 
form- factor calculations [H, @, 0,0] to determine \V u b\. 

The measurement uses a sample of approximately 232 
million BB pairs, corresponding to an integrated lumi- 
nosity of 211 fb — 1 , recorded near the T(4S) resonance 
with the BABAR detector Q at the PEP-II asymmetric- 
energy e + e~ storage rings. We use a detailed Monte 
Carlo (MC) simulation to estimate the signal efficiency 
and the signal and background distributions. 

In the first analysis, we reconstruct -Btag in the semilep- 
tonic decay B — > D^*>lv. We reconstruct D° mesons 
in K~ir + , K~it + it + 'n;~ , _ftT~7r + 7r°, and K®ir + ir~ decays, 
and D + mesons in K~Tr + ir + decays 10]. The D mass 
resolution (a) is between 4.6 and 12.9 MeV depending 
on the decay channel. The mass of the D candidate 



is required to be within 2.6a and 3.0c of the expected 
value for the B° and B + channels, respectively. We also 
use a sideband sample, in which the D candidate mass 
is more than 3ct away from the nominal value, for sub- 
tracting the combinatoric background. We reconstruct 
D* + mesons in D°ir + and D + ir° decays. The mass dif- 
ference between the D* and D is required to be within 
3 MeV of the expected value [ll| . The reconstructed D 
and D* candidates are paired with a charged lepton with 
a center-of-mass (cm.) momentum |p^| > 0.8 GeV to 
form a Y = D^*H system. If the D decay contains a 
charged kaon, the lepton must have the same charge as 
the kaon. The lepton and the D meson are required to 
originate from a common vertex. Assuming that only 
a massless neutrino escaped detection, we calculate the 
cosine of the angle between the B and Y momenta as 
cos9 BY = (2E B Ey -m 2 B -m%)/(2\p B \\VY\)> where m B , 
my, E B , Ey, Pb, Py refer to the masses, cm. energies, 
and momenta of the B and Y, respectively. For back- 
ground events, cos by does not correspond to the cosine 
of a physical angle and can extend outside ±1. We apply 
a loose selection of | cos B y | < 5 at this stage. 

After identifying the -Btag meson, we require the re- 
maining particles in the event to be consistent with a 
B — » -nlv decay. Charged tracks that are not identi- 
fied as a lepton or a kaon are considered charged pion 
candidates. Neutral pion candidates are formed from 
pairs of photon candidates with invariant mass between 
115 and 150 MeV. For the B° channel, the lepton must 
have Ip^l > 0.8 GeV, and its charge must be oppo- 
site to that of the charged pion. The lepton charge 
must be opposite to that of the -B tag for the B + chan- 
nel. We reject the lepton candidate if, when combined 
with an oppositely-charged track, it is consistent with a 
J/tp — > £ + £~ decay or a photon conversion. Once the sig- 
nal B candidate is identified, we require that the event 
contain no other charged particles and small total cm. 
energy E ICS of the residual neutral particles. In mea- 
suring E res , we remove the neutral candidates that are 
consistent with coming from a D* — * Dir° or D'y de- 
cay, bremsstrahlung from an electron, or beam-related 
background. We require B rcs < 70 MeV for the B° chan- 
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nel and E les < 250 MeV for the B + channel, the latter 
being relaxed to allow for additional photons from de- 
cays of D*° and higher resonances. We calculate the 
cosine of the angle between the B and ir£ momenta as 
cos6 B7T i = (2E B E^i - m% - m^)/(2|p B ||p^|), where 
m^e, E-xi, Ptti are the mass, cm. energy, and momentum 
of the n£ system, respectively. We require | cos^s^l < 5. 

Ignoring the small cm. momentum of the B meson, 
the invariant mass squared of the lepton-neutrino system 
in a B — ► tt£v decay can be inferred as q 2 — (tub — 
E^) 2 — |pvr | 2 , where E n and p,r are the cm. energy and 
momentum of the pion. We divide the data into three 
bins: q 2 < 8GeV 2 , 8 < q 2 < 16GeV 2 , and q 2 > 16GeV 2 . 
We use simulated B ir£v events to estimate and to 
correct for the small (< 8%) migration between the q 2 
bins due to resolution, which is approximately 0.8 GeV 2 
at q 2 = 8 GeV 2 and improves with increasing q 2 . 

Having identified the two B mesons that decayed 
semileptonically, conservation of the total momentum de- 
termines the angle <f>B between the direction of the B mo- 
menta and the plane defined by the Y and ir£ momenta: 

2 , COS 2 9 BY + COS 2 6 B-Rg + 2 COS 8 BY COS 9 Btt£ COS 
COS <pB = o 

sin 7 

(1) 

where 7 is the angle between the Y and ir£ momenta. 
The variable cos 2 4>b satisfies cos 2 0s < 1 for correctly 
reconstructed signal events, and is broadly distributed 
for the background (see Fig. [TJ). We use the cos 2 4>b dis- 
tributions to extract the signal yield in the data in each 
q 2 bin. We did not require stringent cuts on cos 9by and 
cos 9 bttI because they are incorporated in cos 2 (f> B - 

We express the data distribution as a sum of three 
contributions: dN/ d cos 2 4>b = N s i g V s i g + NbkgPbkg + 
N cm bP cm b, where N c and V c are the number of events 
and the probability density function (PDF) for each cat- 
egory c, defined as the signal (sig), background with 
correctly- reconstructed D mesons (bkg) , and other back- 
grounds (cmb). The events in the D mass sideband are 
also used in the fit to constrain the -/V cm b'Pcmb term. 
The PDF shapes are determined from the MC simula- 
tion. The signal PDF is a combination of a smeared 
step function and an exponential tail. The background 
PDFs are either an exponential plus constant or a sec- 
ond order polynomial. The two data samples (D mass 
peak and sideband) and the MC samples are used in an 
unbinncd maximum likelihood fit that determines iV s i g , 
-^bkg) -^Vcmb, and the PDF parameters simultaneously. 
Figure [T] shows the fit results summed over the q 2 bins. 
We find the signal yields and their statistical errors to 
be 57±i2 events and 92± 2 ,| events for the B° and B+ 
channels, respectively. 

We use simulated B — > -kIv events to estimate the 
signal efficiencies. Control samples are used to derive 
corrections for the data-MC differences in the B tag re- 
construction, charged and neutral particle reconstruc- 




cos 2 (|) b cos 2 ([» b 



FIG. 1: Distributions of cos 2 </>b of the a) B° — » tt~£ + i> and b) 
B + — » tv°£ + u candidates. The points with error bars and the 
shaded histograms are the data in the D mass peak and side- 
band, respectively. The curves are the fit results representing 
the total (solid), background (dashed), and 'cmb' (dotted) 
components defined in the text. The fits were performed in 
bins of q 2 , but the results shown are for the complete q 2 range. 



tion, and lepton identification. The largest uncertainty 
comes from the Btag reconstruction efficiency, which is 
determined from a sample of events in which two non- 
overlapping Stag candidates are reconstructed. The effi- 
ciency correction factors for the Stag reconstruction are 
found to be 1.00 ± 0.07 and 0.99 ± 0.02 for the B° and 
B + channels, respectively. The average signal efficiencies 
after the correction are 1.1 x 10~ 3 for the B° channel 
and 3.0 x 10~ 3 for the B + channel. The latter is larger 
mainly because of the higher efficiency of reconstructing 
a D° meson compared with a D + or D* + meson. 

The measured branching fractions are summarized in 
Table |H The largest sources of systematic error [l2j are: 
the -Btag reconstruction efficiency (discussed above), the 
shape of the background cos 2 4>b distribution (studied 
with control samples that fail the signal selection crite- 
ria), and the branching fractions of the B semileptonic 
decays other than B — > tt£v (varied within the current 
knowledge [llj). 

In the second analysis, we reconstruct the -Btag meson 
in a set of purely hadronic final states B — > D^*>X. We 
reconstruct D° mesons in-ftT - 7r + , K~ir + ir°, K~Tr + ir + ir~ , 
and K®Tr + ir~ decays, and D + mesons in K~7r + -K + , 
A'-tt+tt+tt , K°tt+, K°7T+TT°, and K° s tt+ tt+ ^ decays. 
The D* mesons are reconstructed in D°ir + , D°ir°, and 
_D°7 decays. The hadronic system X has a total charge 
±1 and is composed of riiir^ +n2K ± +71311° +ri4K® where 
ni +ri2 < 6, n.3 < 3 and 77,4 < 3. The total reconstruction 
efficiency for a B° (B+) meson is 0.3% (0.5%). 

We separate correctly-reconstructed B tag mesons from 
the background using two kinematic variables: the beam- 
energy substituted mass tties = V / s/4 — |ps| 2 and the 
energy difference AE = Eb — y/s/2, where y/s is the cm. 
energy of the e + e~ system. We select signal candidates 
in mode-dependent AE windows around zero. We apply 
a loose selection 5.2 < m^s < 5.3 GeV and fit the tties 
distribution at a later stage to extract the signal yield. 
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TABLE I: Partial and total branching fractions, in units of 10 4 , measured with the semileptonic and hadronic tag analyses. 
The q ranges are in GeV 2 . The errors are statistical and systematic. The combined results are expressed as B° — > it £ + v 
branching fractions. 

q 2 < 8 8 < q 2 < 16 q 2 > 16 q 2 < 16 Total 

B° Semileptonic 0.50 ± 0.16 ± 0.05 0.33 ± 0.14 ± 0.04 0.29 ± 0.15 ± 0.04 0.83 ± 0.22 ± 0.08 1.12 ± 0.25 ± 0.10 
Hadronic 0.09 ± 0.10 ± 0.02 0.33 ± 0.15 ± 0.05 0.65 ± 0.20 ± 0.13 0.42 ± 0.18 ± 0.05 1.07 ± 0.27 ± 0.15 
Average 0.38 ± 0.12 ± 0.04 0.33 ± 0.10 ± 0.03 0.47 ± 0.13 ± 0.06 0.72 ± 0.16 ± 0.06 1.19 ± 0.20 ± 0.10 
B + Semileptonic 0.18 ± 0.08 ± 0.02 0.45 ± 0.13 ± 0.05 0.10 ± 0.12 ± 0.04 0.63 ± 0.16 ± 0.06 0.73 ± 0.18 ± 0.08 
Hadronic 0.16 ± 0.11 ± 0.03 0.39 ± 0.16 ± 0.06 0.26 ± 0.12 ± 0.06 0.56 ± 0.19 ± 0.08 0.82 ± 0.22 ± 0.11 
Average 0.18 ± 0.07 ± 0.02 0.43 ± 0.10 ± 0.04 0.22 ± 0.09 ± 0.05 0.61 ± 0.12 ± 0.05 0.82 ± 0.15 ± 0.09 
Combined 0.36 ± 0.09 ± 0.03 0.52 ± 0.10 ± 0.04 0.46 ± 0.10 ± 0.06 0.87 ± 0.13 ± 0.06 1.33 ± 0.17 ± 0.11 



After reconstructing the Stag) we look for the signa- 
ture of a B — > tt£v decay in the recoiling system. The 
selection criteria for the pion and lepton candidates are 
similar to the first analysis, except a) the minimum |p^| 
for electrons is 0.5 GeV, and b) the ir° mass window is 
110-160 MeV. We require E ICS < 450 MeV for the B° 
channel to reduce the B° — > p~t + v background, and no 
requirement is made for the B + channel. 

The full reconstruction of -Btag allows us to determine 
the neutrino four-momentum precisely from the missing 
four-momentum p miss = p T (4S) - PB tas - P-k - Pi- The 
missing mass squared w^iss peaks near zero for the signal 
and extends above zero for the background (see Fig. ^ . 
We require Iw^issl < 0.3 GeV 2 for the B° channel and 
-0.5 < TO 2 niss < 0.7 GeV 2 for the B + channel, with the 
latter being broader and asymmetric due to the resolu- 
tion of the 7T° energy measurement. 

Precise knowledge of p m i ss allows us to calculate q 2 
with small uncertainties. We divide the signal candi- 
dates into the same three q 2 bins as before, and subtract 
the small bin-to-bin migration as background. In each 
q 2 bin, we obtain the number of correctly-tagged events 
by an unbinned maximum likelihood fit to the tties dis- 
tribution. The PDF for the signal is determined from 
MC simulation as a Gaussian function joined to an ex- 
ponential tail. For the background, we use a threshold 
function of the form x\/l — x 2 exp(— £(1 — x 2 )), where 
x = 2toes/v / * an d the parameter £ is allowed to float in 
the fit. Fig. [2] shows the distribution obtained by 

splitting the data samples in bins of and repeating 

the 777.ES fit. 

The signal side of the correctly-tagged events may 
not be a B — ► nlv decay. Contributions from this 
type of background are estimated with the MC simula- 
tion, as indicated by shaded histograms in Fig. [51 which 
are scaled to match the data in the sideband region 
1 < 777 2 niss < 4 GeV 2 . After background subtraction, wc 
find signal yields of 31 ± 7 events and 26 ± 7 events for 
the B° and B + channels, respectively, where the errors 
are statistical. 

Instead of estimating the absolute signal efficiency, we 




FIG. 2: Distributions of m^ iss of the a) B° — > tt~£ + v and 
b) B + — > tt°£ + v candidates. The points with error bars are 
the data. The histograms represent, from the lightest to the 
darkest, the MC simulation of the B — > tt£v signal, b — > uiv, 
b — » civ, and other backgrounds. The arrows indicate the 
regions in which the signals are extracted. 

normalize the signal yield to the number of inclusive B 
semileptonic decays, B — > Xlv, in the recoil of Stag- The 
reconstruction efficiencies of the -Btag and of the lepton 
cancel to first order in the ratio between the yields of the 
signal and normalization samples. The inclusive branch- 
ing fraction B(B -> Xlv) is taken as 10.73 ± 0.28% \}^. 
The yield of the normalization sample is extracted by a 
fit to the 777es distribution. The component of the back- 
ground that peaks in the tties distribution is estimated 
from the MC simulation and subtracted. Efficiency dif- 
ferences between the signal and normalization samples 
are estimated with the MC simulation, and the corre- 
sponding corrections are applied to the result. 

The measured branching fractions are summarized in 
Tabled] The largest source of systematic error is the lim- 
ited statistics of the signal MC sample. Other significant 
sources include the modeling of the signal PDF (studied 
with alternative fitting methods), photon-energy mea- 
surement, 7T° reconstruction, muon identification, and 
the branching fractions of non-signal B — > X u £v decays. 

We take weighted averages of the measured partial 
branching fractions in each q 2 bin. The results for the B° 
and B + channels are consistent with the isospin relation 
T(B° -> ir-£+v) = 2T(B+ -> ir°£+v) and the lifetime 
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TABLE II: Values of \V u b\ derived using the form factor cal- 
culations. The first two errors on \V u b\ come from the sta- 
tistical and systematic uncertainties of the partial branching 
fractions. The third errors correspond to the uncertainties on 
AC due to the form-factor calculations, and are taken from 
Refs.lSQ! 





(GeV 2 ) 


AC 


(PS" 1 ) 




\v ub \ (io- 3 


) 


Ball-Zwicky [5] 


< 


16 


5.44 


± 1.43 


3.2 


±0.2 ±0.1 


+ 0.5 
-0.4 


Gulez et al. [6] 


> 


16 


1.46 


±0.35 


4.5 


±0.5 ±0.3 


+0.7 
-0.5 


Okamoto et al. [7] 


> 


16 


1.83 


±0.50 


4.0 


±0.5 ±0.3 


+0.7 
-0.5 


Abada et al. [8J 


> 


16 


1.80 


±0.86 


4.1 


±0.5 ±0.3 


+ 1.6 
-0.7 



ratio t b +/t b o = 1.081 ± 0.015 [ll|, with X 2 = 5.2 for 
3 degrees of freedom. Assuming isospin symmetry, we 
combine the B° and B + channels and express the results 
as the B° branching fraction in the last row of Table Q] 
The overall x 2 is 10.2 for 9 degrees of freedom. 

We extract \V u b\ from the partial branching fractions 
AB using JK b | = y/AB/(T B oA£), where t b o = (1.536 ± 
0.014) ps [lH is the B" lifetime and AC = Ar/|Kb| 2 is 



the normalized partial decay rate predicted by the form- 
factor calculations. We use the light-cone sum rules cal- 
culation 5(1 for q 2 < 16 GeV 2 and the lattice QCD calcu- 
lations [aSH for I 2 > 16 GeV 2 . The results are shown 
in Table OH 

In conclusion, we have measured the B — > ir£v branch- 
ing fraction as a function of q 2 using tagged B meson 
samples, and have extracted \V u b\. The measured total 
branching fraction, B(B° -> 7T-£ + v) = (1.33 ± 0.17 sta t ± 
0.11 S y S t) x 10~ 4 , has the smallest systematic uncertainty 
among the existing measurements 0, [3, 0] thanks to the 
superior signal purity, and the overall precision is com- 
parable to the best. Using theoretical calculations of the 
form factor, we obtain values of \V u b\ ranging between 
3.2 x 10~ 3 and 4.5 x 10~ 3 . As an example, the recently 
published unquenched lattice QCD calculation [§[ gives 
|Kfe| = (4.5 ± 0.5 stat ± 0.3 syst ±°;£ PF ) x 10" 3 . Improve- 
ment will be possible with additional data combined with 
more precise form-factor calculations. 
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TABLE I: Fractional systematic errors (in %) of the measured partial branching fractions. The q 2 bins are defined as 1: 
q 2 < 8GeV 2 , 2: 8 < q 2 < 16GeV 2 , and 3: q 2 > 16GeV 2 . The + and x symbols indicate if the error is additive (+) or 
multiplicative (x). 





Stag 


B° 


semilep. 


B 4 


semilep. 


B° 


hadronic 


B+ 


hadronic 




bin 


1 


2 


3 


1 


2 


3 


1 


2 


3 


1 


2 


3 


B —> nlv form factor 


X 


1.0 


0.5 


1.1 


0.9 


0.5 


4.5 


0.3 


0.2 


0.1 


0.3 


0.2 


2.2 


B — > XJLv background 


+ 


1.9 


2.9 


3.8 


2.0 


3.5 


7.7 


0.2 


0.2 


0.2 


2.6 


2.6 


2.6 


B — » X u lv background 


+ 


0.8 


1.7 


6.9 


1.2 


1.7 


12.1 


4.2 


4.2 


4.2 


1.7 


1.7 


1.7 


B(B -> Xiv) 


x 




not 


applicable 




2.6 


2.6 


2.6 


2.6 


2.6 


2.6 


B(Y(4S) -> B°B°) 


X 


1.6 


1.6 


1.6 


1.6 


1.6 


1.6 




not applicable 




Final-state radiation 


X 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


1.2 


Stag efficiency 


X 


7.3 


7.3 


7.3 


4.3 


2.5 


12.9 


0.7 


0.7 


0.7 


1.4 


1.4 


1.4 


q 2 resolution 


X 


1.6 


1.3 


1.2 


1.2 


4.5 


18.0 






negli 


gible 






Fit method 


+ 


1.4 


2.1 


5.7 


4.8 


6.2 


32.8 


5.7 


5.7 


5.7 


2.7 


2.7 


2.7 


Lepton identification 


X 


1.6 


1.9 


1.9 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


Charged track reconstruction x 


1.6 


1.6 


1.6 


0.8 


0.8 


0.8 


1.1 


1.1 


1.1 


1.4 


1.4 


1.4 


Neutral energy reconstruction x 


negligible 


3.2 


3.0 


6.3 


1.2 


1.2 


1.2 


3.7 


3.7 


3.7 


Number of BB events 


X 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 




not applicable 




MC statistics 


X 


5.2 


5.1 


4.6 


4.7 


4.2 


6.8 


18.3 


11.8 


17.6 


19.8 


14.7 


23.0 


Total 


10.0 


10.4 13.6 


9.7 


10.9 


43.5 


20.1 


14.4 


19.4 


21.0 


16.3 


24.1 



TABLE II: Values and errors (in unit of 10 -4 ) of the combined partial branching fractions. The errors are separated into statis- 
tical, multiplicative systematic, and non-multiplicative systematic components, and the covariance matrices for the systematic 
components are given. The q 2 bins are defined as 1: q 2 < 8GeV 2 , 2: 8 < q 2 < 16GeV 2 , and 3: q 2 > 16GeV 2 . 



bin 


12 3 


Value 


0.355 0.518 0.457 


Statistical error 


0.086 0.097 0.104 


Multiplicative systematic error 


0.028 0.038 0.052 


Covariance 1 

2 
3 


1.000 0.593 0.471 
0.593 1.000 0.504 
0.471 0.504 1.000 


Non-multiplicative systematic error 


0.008 0.017 0.021 


Covariance 1 

2 
3 


1.000 0.986 0.791 
0.986 1.000 0.881 
0.791 0.881 1.000 



